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Summary 

The near-field noise characteristics of a propfan 
operating at an angle of attack are examined utilis- 
ing the unsteady pressure field obtained from a three- 
dimensional Euler simulation of the propfan flowfield. 
The near-field noise is calculated employing three dif- 
ferent procedures: a direct computation method in 
which the noise field is extracted directly from the 
Euler solution, and two acoustic- analogy-based fre- 
quency domain methods which utilize the computed 
unsteady pressure distribution on the propfan blades 
as the source term. The inflow angles considered are 
-0.4, 1.6, and 4.6 degrees. The results of the direct 
computation method and one of the frequency do- 
main methods show qualitative agreement with mea- 
surements. They show that an increase in the inflow 
angle is accompanied by an increase in the sound pres- 
sure level at the outboard wing boom locations and a 
decrease in the sound pressure level at the (inboard) 
fuselage locations. The trends in the computed az- 
imuthal directivities of the noise field also conform to 
the measured and expected results. 

Introduction 

The advanced propfan program initiated by 
NASA in the mid-seventies culminated in flight 
test /demonstration of single and counter-rotation 
propfan systemsfl]. The NASA Lewis advanced tur- 
boprop project included wind tunnel tests of a series 
of propfan models, and wind tunnel and flight tests 
of a large-scale propfan system. A large-scale 9-foot 
(2.74m) diameter, NASA-Hamilton Standard single 
rotation eight bladed tractor configuration was flight 
tested in 1987-88. The propfan was mounted on the 
left wing of an instrumented Gulfstream II testbed 
aircraft, modified for the propfan tests (Fig. la). The 
Propfan Test Assessment program (PTA) was con- 
ducted by Lockheed Aeronautical Systems Company 
for NASA Lewis Research Center. The overall 


objectives of the PTA program were to evaluate the 
propfan structural integrity, source noise, cabin noise, 
and enroute noise[2,3]. 

In the PTA program, a nacelle tilt arrangement 
was employed to vary the inflow angle to the prop- 
fan. The nacelle tilt angles considered were -3, -1 
(tilt down) and +2 (tilt up) degrees. The average in- 
flow angle to the propfan, however, is dependent not 
only on the nacelle tilt angle but also on the airplane 
angle of attack and the upwash angle at the propfan 
(Fig. lb). For the cruise operating conditions of in- 
terest here the airplane mean angle of attack was 1.6 
degrees. The upwash angle at the propfan estimated 
from panel method calculations^] was about 1.0 de- 
gree. Thus, the estimated propfan inflow angles[4] 
corresponding to the nacelle tilt angles of -3, -1 and 
+2 degrees were -0.4, 1.6 and 4.6 degrees, respectively. 

The variations in inflow conditions were used to 
evaluate the effects of cyclic stresses on the propfan 
performance over a wide range of operating condi- 
tions. In particular, the effects of inflow conditions 
on the near-field noise characteristics were studied for 
takeoff and cruise, operating conditions. In this paper 
attention is restricted to the near- filed noise charac- 
teristics at cruise conditions. The measurements show 
that an increase in inflow angle is accompanied by 
an increase in sound pressure level (SPL) at the wing 
boom and a simultaneous decrease in SPL at the fuse- 
lage^]. In the PTA tests and predictions reported in 
this paper, the sense of rotation of the propfan is in- 
board up. 

The prediction of the effect of inflow angle on 
the near-field SPL characteristics requires knowledge 
of the unsteady flowfield of the propfan. The numeri- 
cal solution of the three-dimensional Euler equations 
governing the flowfield offers an enticing possibility. 
Nallasamy and Groeneweg[5] reported unsteady Eu- 


ler solutions of the flowfield of a propfan at angle of 
attack for three inflow angles of -0.4 ,1.6 and 4.6 de- 
grees. The conditions chosen for that study corre- 
spond to the PTA flight test cruise conditions with a 
blade tangential tip speed of 797 ft/sec (243 m/sec). 
The computed unsteady blade surface pressure dis- 
tributions were compared with wind tunnel data[6,7] 
and found to be in reasonable agreement. 

In the present investigation, the near-field noise 
characteristics are examined as a function of inflow 
angle employing the unsteady flowfield solutions of 
Nallasamy and Groeneweg[5]. The near-field SPLs 
are computed using three different methods to be de- 
scribed later. The computed SPLs are compared with 
the PTA measurements. 

Unsteady flowfield solutions 

The unsteady three-dimensional Euler solutions 
reported in reference 5 are used for the present near- 
field acoustic calculations. The salient features of the 
unsteady solutions are briefly described here for clar- 
ity and aid in understanding of the acoustic calcula- 
tions. Fig. 2 shows the computed azimuthal varia- 
tions of blade loading due to angular inflow for three 
inflow angles, -0.4 , 1.6 and 4.6 degrees. The azimuth 
angle is measured from the vertical(top-dead-center, 
Fig. lc). The sinusoidal variation of the integrated 
blade loading due to angular inflow is evident. A 
substantial variation of the loading(±81%) about the 
mean is observed for 4.6 degrees inflow angle. The 
variations of the integrated blade loading with inflow 
angle is shown in Fig. 3 for two azimuthal locations 
<f> = 90 degrees (boom side), and <f> = 270 degrees 
(fuselageside) . An increase in the inflow angle causes 
the blade loading on the boom side to increase, but 
that on the fuselage side to decrease. 

The variation of elemental power coefficient with 
inflow angle for four azimuthal positions, <f> = 0, 90, 
180, and 270 degrees is shown in Fig. 4. At <j> = 0 
and 270 degrees locations the blade loading decreases 
with increase in the inflow angle throughout the blade 
span. The blade loading increases with increase in the 
inflow angle over the entire span at <f> = 90 and 180 
degrees locations. 

Typical azimuthal variations (waveforms) of un- 
steady pressure coefficients are shown in Fig. 5. The 
figure shows the waveforms at the 75% radial station 
at three normalized chordwise locations, 0.1, 0.5 and 
0.9, on the suction and pressure sides of the blade. 
The blade response increases with inflow angle as 
would be expected. Non-sinusoidal responses are ob- 
served for 4.6 degrees inflow angle both on the suc- 
tion and pressure sides due to the presence of shock 


waves[5]. The variation of blade response with ax- 
ial distance (chord) is also illustrated. These blade 
surface pressure distributions form the input for the 
two frequency-domain acoustic calculations described 
later. Further details of the azimuthal, chordwise, 
and spanwise variations of the blade surface pressures 
and blade loading may be found in [4]. 

Near-field Acoustic Calculations 

The noise levels at the microphone locations on 
the boom and fuselage sidelines, at a distance of 1.12D 
(D is the propfan diameter) from the axis, were com- 
puted for comparisons with measured values. The 
azimuthal and axial locations of the boom micro- 
phones and fuselage microphones opposite the boom 
are shown in Figs. 6a and 6b, respectively. Note that 
the geometric observation angle at any microphone 
changes with the nacelle tilt since the propfan cen- 
ter moves up or down relative to the fixed boom and 
fuselage microphones(Fig. 6a). 

The noise predictions were obtained employing 
three methods: a direct numerical computational pro- 
cedure and two frequency domain analysis methods. 
A brief description of each of these methods is given 
below. 

(i) Direct numerical computation procedure (DCP) 

In this method the propfan near-field noise is ex- 
tracted from the pressure field obtained from solving 
the three-dimensional Euler equations. For a steady 
flow(Le., zero inflow angle) since the flowfield in each 
blade passage is identical, the Euler solution is usu- 
ally carried out for only one blade passage. From 
such an Euler solution, a spatial pressure distribution 
in the azimuthal (blade to blade) direction is obtained 
for a specified observer location. The spatial pressure 
pattern is transformed into the time domain using 
the rotational speed and blade count of the propeller. 
Finally, the pressure-time history thus obtained is 
Fourier analyzed for its spectral content. 

Since in the DCP method the acoustic pressure is 
extracted directly from the field pressure, it includes 
the nonlinear propagation effects which are not typi- 
cally accounted for in conventional methods based on 
acoustic analogy. The DCP method, however, suf- 
fers from the limitations of the flowfield resolution 
in the regions away from the propfan blade. This 
procedure (DCP) was previously employed for steady 
flow (i.e., zero inflow angle) conditions by Korkan et 
al.[8j and Whitfield et al.[9]. Their Euler solution 
suffered from low grid resolution in regions away from 
the blade due to grid clustering schemes typically em- 
ployed in these computations. Hence, Whitfield et al. 
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computed the sound pressure levels using DCP at dis- 
tances close to the blade and then extrapolated their 
results to the observer (microphone) locations of in- 
terest. 

In contrast to the zero inflow angle case described 
above, for a propeller operating at an angle of attack, 
the unsteady Euler equations must be solved for the 
entire flowfield (i.e., all blade passages). The flow- 
field in each passage is different due to variations in 
local flow incidence angles as a result of the angular 
inflow. Also, in this case the pressure-time history at 
any observer position can be obtained directly from 
the unsteady three-dimensional Euler solution. The 
pressure waveform is obtained by storing the pressure 
computed at an observer location as the blades rotate 
past it as is done in a physical experiment. But no 
ensemble averaging is done in the computations! The 
harmonic content is then extracted as in the steady 
case by a Fourier representation of the pressure wave- 
form. 

In the unsteady Euler computations of PTA 
cruise test cases for the three inflow angles, a grid 
distribution of the type shown in Fig. 7 was em- 
ployed. This figure shows the axial-radial surface grid. 
The waveforms extracted at the PTA sideline micro- 
phone location(1.12D) indicated that the flow resolu- 
tion there is not sufficient for such purposes. Hence, 
following the procedure of Whitfield et al.[9], as a 
first approximation, the acoustic pressure waveforms 
are computed at a radial location of 0.53D from the 
propfan axis. The observer angles/axial locations of 
the PTA microphones and the locations selected for 
computation are shown in Fig. 7. The computed re- 
sults are then extrapolated to the actual PTA micro- 
phone sideline distance(1.12D) using a simple inverse 
square distance law. 

Unlike the DCP method described above, the fre- 
quency domain methods utilize the concept of acous- 
tic analogy to estimate the acoustic radiation from 
rotating propellers. In these methods the propfan 
blade unsteady surface pressure distributions are used 
as source strengths which then replace the propfan al- 
together. 

(ii) Hansons’s Frequency- Domain Theory (FDH) 

This theory known as the helicoid al surface the- 
ory, was developed by Hanson [10]. It is widely used 
for propfan noise studies. It employs the helicoidal 
surface representation of the moving sources to derive 
a near/far field estimate of the noise radiation from 
propfans. Based on this theory and its companion 
aerodynamic theory [11] a computer program called 


Unified Aero-Acoustic Program (UAAP) was devel- 
oped by Hamilton Standard for NASA Lewis. In the 
present study the acoustics part of UAAP was used 
to compute the near-field noise for the PTA angular 
inflow cases. The azimuthal variation of blade loading 
(waveform) obtained from the Euler solution is rep- 
resented by a Fourier series and the unsteady loading 
harmonics are used as input to the UAAP. The UAAP 
computes the SPLs, taking one blade loading order 
(harmonic) input at a time. Once noise levels from 
all loading orders of interest are computed, they are 
added together accounting for their proper phases. 
The SPLs reported here are the ones computed for 
the first three unsteady blade loading orders. 

(iii) frequency Domain Formulation of Envia (FDE) 

This method uses a version of the Ffowcs 
Williams-Hawkings[12] equation derived for a uni- 
formly moving medium. This formulation is a natural 
choice for comparisons with the PTA inflight mea- 
surements as well as wind tunnel data. Like the FDH 
method, FDE method is carried out in the frequency 
domain but unlike the former, it does not rely on the 
near/far approximations of a certain integral. In the 
FDE method this integral is evaluated utilizing an 
efficient numerical integration scheme (Romberg in- 
tegration). The use of a simple numerical integra- 
tion scheme further affords the flexibility of using a 
source strength that can have an arbitrarily complex 
time-dependent behavior without the need to decom- 
pose the source strength time history into its con- 
stituent Fourier components. Furthermore, the use 
of the numerical integration allows the inclusion of 
the cross-flow convective phase effects (i.e., Mani’s 
“added phases* [13] ) which are present whenever the 
inflow angle into the propeller is non-zero. In fact, 
in the FDE method, this effect is accounted for in an 
exact manner for arbitrarily large cross-flow compo- 
nents. The contribution due to Mani’s added phases 
to the computed SPL is, however, small for the in- 
flow angles and cruise operating conditions considered 
here. 

In addition to the above issue, there is another 
important difference between the FDH and FDE 
methods. In the FDH method, the local integrated 
“pressure force* (i.e., lift) is taken to be perpendic- 
ular to the local direction of motion rather than the 
local blade surface. This “2-D” description of loading 
impacts the acoustic predictions in two ways. Firstly, 
since it basically excludes any radial loading contri- 
bution it doe s not account for the effect of the radial 
dipole which cannot be neglected in the vicinity of 
the plane of rotation of the propeller. Secondly, by 


3 



placing the sectional lift perpendicular to the local 
direction of motion, all of the local acoustic dipoles 
at a given radial station are, in effect, aligned with 
each other. This, in turn, tends to reduce the mutual 
interference between them which otherwise might ex- 
ist. By placing the local loading perpendicular to the 
local blade surface, as is done in the FDE method, 
a more realistic distribution for the acoustic sources 
is obtained. It should be mentioned that the effects 
alluded to here tend to be unimportant for the con- 
ventional straight blade propellers of the type used in 
general aviation. Details of the FDE formulation may 
be found in [14]. 

Results and Discussion 

The near-field noise levels of the propfan oper- 
ating at angle of attack have been computed for the 
PTA cruise test conditions. The propfan inflow angles 
considered are -0.4, 1.6, and 4.6 degrees. The cruise 
test conditions are: Mach number — 0.8 , advance ra- 
tio = 3.12, altitude = 35000 ft (10,670 m) and blade 
tangential tip speed = 797 ft/sec (243 m/sec ). The 
near-field noise levels computed by the three methods 
described above are compared with the PTA measure- 
ments. 

The wing acoustic boom on the PTA tapers from 
6-1/4 in. (15.8 cm) diameter at the microphone at 
1.0D (closest to the wing) to 1-5/8 in. (4.1cm) di- 
ameter at -0.5D (forward) microphone. The scatter- 
ing from the boom affects the SPLs measured at the 
boom microphones. The boom scattering at the mi- 
crophone locations computed using reference 15 was 
reported by Dittmar[l6]. The pressure amplifications 
for the boom microphones going from forward to aft 
are 0, 0.5, 0.5, 1.0 and 1.5 dB for the blade passing 
frequency (BPF) tone; 0*5, 1.5, 1.5, 2.0, and 2.5 dB 
for 2BPF; and 1.0, 2.0, 2.5, 3.5, and 4 dB for 3BPF. 
These boom amplifications are subtracted from the 
measured boom tones (Figs. 8,10 and 11) for com- 
parisons with the computed results. The fuselage is 
assumed to have a 6dB pressure amplification and this 
value is subtracted from the fuselage data (Fig. 9) for 
comparison with predictions. 

Figure 8 shows the effect of inflow angle on the 
area-maximum SPL at the blade passing frequency at 
the outboard wing boom microphone side-line. The 
area-maximum SPL here is the maximum of the SPLs 
recorded at five boom microphone locations (see Fig. 
6b for microphone locations). At these locations, 
measurements (Fig. 8) s how that an increase in in- 
flow angle increases the fundamental tone SPL. The 
SPLs computed by the DCP and FDE methods show 


the same trend. On the boom side, the SPL increases 
with increase in inflow angle as does the blade loading 
(Fig. 3). The results of the FDH method show little 
sensitivity to the inflow angle. Such an insensitivity 
in the predictions from the FDH method was also re- 
ported in[l7] (see also[4j). It should be noted that the 
predictions reported in[l7] were obtained using a dif- 
ferent unsteady blade loading distribution. In view 
of the discussion in the previous section concerning 
the differences between the FDH and FDE methods, 
perhaps the apparent insensitivity of the FDH pre- 
dictions to the inflow angle changes is related to the 
manner in which the loading dipoles are implemented 
in that method. 

Figure 9 shows the effect of inflow angle on the 
area-maximum SPL at the BPF at the inboard fuse- 
lage microphone sideline. The area-maximum SPL 
here refers to the maximum of the SPLs recorded 
at the eight fuselage microphone locations (Fig. 6b 
shows these microphone locations). On the fuselage 
side (Fig. 9), the measurements show that an increase 
in inflow angle causes a decrease in fundamental tone 
SPL. The results of DCP and FDE calculations show 
a similar trend. However, in the computations the 
slope of the line is steeper than in the measurements. 
On the fuselage side, the SPL decreases with increases 
in inflow angle as does the blade loading (see Fig. 3). 
Again, the results from the FDH method show lit- 
tle sensitivity to inflow angle changes. It should be 
noted that the computed results shown in Fig. 9 are 
free field and do not include corrections for bound- 
ary layer refraction and wing reflection effects. Since 
no model for these effects has been included in the 
present computations further comparisons with data 
will be limited to boom microphone locations. 

The predicted sensitivities of the first, second 
and third harmonic SPLs to inflow angle by the three 
methods at the boom microphone locations are com- 
pared with measurements in Figs. 10a - lOe. The 
fundamental tones measured at all the five boom mi- 
crophone locations generally show an increase with 
inflow angle. This trend is well represented by the two 
predictions, DCP and FDE. The computed absolute 
magnitudes, however, show a mixture of under- and 
over- predict ion. The DCP computations depend on 
the flow field resolution at the computational points. 
Since the location of the forward microphone at -0.5D 
is in the region of poor flow field resolution, no results 
of DCP for this location are presented here and in Fig. 
11. The FDH calculations show a marginal increase 
of SPL with i nflo w angl e at t he two microphones lo- 
cated far aft (l.OD) and forward (-0.5D). Th^ mea- 
8urements may have been influenced by scattering at 
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the boom surface and reflections from the wing both 
of which depend on the microphone location. The 
measurements have been corrected only for the boom 
scattering effects as mentioned above. 

The SPL measurements of the second and third 
harmonics at the boom microphone locations do not 
show a consistent trend with inflow angle. The com- 
putations of FDE and DCP show an increase in SPL 
with increase in inflow angle. The results of the DCP 
method, in general, produce levels of higher harmon- 
ics much lower than the measured ones. 

The sideline BPF tone SPL directivities of the 
boom microphones for the three inflow angles, -0.4 
, 1.6 , and 4.6 degrees are shown in Figs. 11a, lib 
and 11c, respectively. In each of the figures, the mea- 
sured maximum SPL occurs at the microphone lo- 
cated 0.25D aft of the plane of rotation (POR). The 
predictions from the DCP and FDE methods also 
show that, in each case, the maximum level occurs 
at the same location. The FDH predictions, however, 
show the maximum level in each case to occur in the 
plane of the propfan. 

Fig. 12 shows the measured and computed sound 
pressure time histories for the boom microphone lo- 
cated at 0.25D aft of the POR for an inflow angle 
of 1.6 degrees. The peak levels computed by FDE 
are higher than those of the DCP predictions and the 
PTA data as implied by the dB levels in Figs. 8-11. 
The computed sound pressure waveforms are similar 
to but differ from the measured one in that the exis- 
tence of a ’zone’ of near-zero pressure amplitude (for 
about 30% of a period) between positive and negative 
peaks in the measurements is not predicted. 

An interesting near-field noise characteristic of 
a propfan operating at angle of attack is the effect 
of inflow angle on the azimuthal directivity. The 
azimuthal directivities computed by FDE and DCP 
methods are shown in Figs, 13a and 13b respectively. 
In both figures, fuselage side (<f> = 270 degrees) and 
above the propfan (<f> = 0 degree) locations show a de- 
crease in fundamental tone SPL with increase in the 
in flow angle. This is consistent with the blade ele- 
mental power coefficient variations in Fig. 4, which 
also show a decrease in blade power with increases 
in inflow angle for <f> = 0 and 270 degrees. On the 
other hand, the boom side ( <f> = 90 degrees) and be- 
low the propfan [<f> = 180 degrees) locations show an 
increase in fundamental tone SPL with increase in 
the inflow angle. This again is consistent with the 
blade power variations in Fig. 4, which show that for 
<l> = 90 and 180 degrees, blade power increases with 


increases in inflow angle. There exist two small az- 
imuthaT regions about <f> = 45 degree and <f> = 225 
degree where the inflow angle variations have little 
effect on the computed fundamental tone SPLs. The 
azimuthal directivity is also consistent with the no- 
tion that with angular inflow (as in takeoff cases) the 
noise level above the propfan decreases and that be- 
low the propfan (on the ground) the noise level in- 
creases! 18, 19]. The minimum noise level predicted by 
the DCP method is much lower than that computed 
with the FDE method. 

Fig. 14 shows the sideline BPF tone directivities 
computed as a continuous function of axial distance 
with the FDE method for boom and fuselage loca- 
tions opposite the boom. On the boom side (Fig. 
14a), an increase in inflow angle increases the funda- 
mental tone SPL over the entire axial distance, -0.5D 
to 1.0D aft of the POR. However, on the fuselage side 
(Fig. 14b), an increase in inflow angle decreases the 
fundamental tone SPL only in the range of axial dis- 
tance 0.0 to 1.0D. At locations ahead of the POR and 
beyond 1.0D aft of it no clear trend is observed. In 
this range, an inflow angle of 4.6 degrees produces 
levels of SPL higher than those for the other two in- 
flow angles. Note that the plot shows only the results 
of the free field computation. 

Concluding Remarks 

The effect of inflow angle on the near-field noise 
of a propfan operating at cruise conditions has been 
studied using a direct computation procedure and 
two frequency domain methods. All three acoustic 
calculations use the same unsteady blade pressure, 
loading and flow field information from an unsteady 
three-dimensional Euler solution. The results of the 
DCP and FDE methods are in qualitative agreement 
with measurements. The computed SPLs at the wing 
boom locations increase with angle of attack whereas 
those on the fuselage locations opposite the boom de- 
crease with increasing inflow angle for an inboard-up 
propfan rotation. Models for the fuselage boundary 
layer refraction, reflections from the wing surface etc. 
need to be included for more meaningful comparisons 
of the computed absolute levels with measurements. 

The direct numerical computation of the noise 
levels from the Euler solution appears to be a useful 
near-field analysis technique. However, the accuracy 
of the direct computation depends on the flowfield 
resolution at the observer position. Advances in com- 
putational fluid dynamics algorithms and availabil- 
ity of supercomputer time and memory are expected 
to make the required field resolution possible. Then, 
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the direct computation procedure may become a vi- 
able and reliable method of computing near-field noise 
characteristics of a propfan. 
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Predicted Unsteady Pressure Coefficient 
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Fig. 5 Computed unsteady blade pressure response at 75% radial station. 
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Fig, 6 PTA aircraft microphone locations selected for predictions. 
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Fig. 7 Computational grid, axial-radial surface and 
points at which SPLs are computed in DCP. 
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Fig. 8 Effect of inflow angle on the area-maximum BPF tone SPL 
at the boom microphone 1.12D sideline. 
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Fig. 9 Effect of inflow angle on the area-maximum BPF tone SPL 
at the fuselage microphone 1.12D sideline, opposite boom. 
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Fig. 10 Effect of inflow angle on the boom microphone SPLs 
for first three harmonics. 
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Fig. 10 Continued. 
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Fig. 10 Concluded. 
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Fig. 12 Computed and measured sound pressure waveforms at 
0.25D aft boom microphone location, inflow angle=1.6°. 



Fig. 13 Computed azimuthal directivity 0.25D aft. 
(a) FDE and (b) DCP predictions. 
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Fig. 14 1.12D sideline directivities at the boom microphone locations 
and fuselage locations opposite boom, FDE predictions. 
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